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Abstract 

A porous material was considered as a platform for optical sensing. It was envisaged that the porous 
material was infiltrated by a fluid which contains an agent to be sensed. Changes in the optical properties 
of the infiltrated porous material provide the basis for detection of the agent to be sensed. Using a 
homogenization approach based on the Bruggeman formalism, wherein the infiltrated porous material 
was regarded as a homogenized composite material, the sensitivity of such a sensor was investigated. For 
the case of an isotropic dielectric porous material of relative permittivity and an isotropic dielectric 
fluid of relative permittivity e*", it was found that the sensitivity was maximized when there was a 
large contrast between and e*"; the maximum sensitivity was achieved at mid- range values of porosity. 
Oh Especially high sensitivities may be achieved for close to unity when e" 2> 1, for example. Furthermore, 

O ■ higher sensitivities may be achieved by incorporating pores which have elongated spheroidal shapes. 

■ 1 Introduction 

X. , , , 

, A simple generic optical sensor may be envisaged as porous material (labelled a, say), vifhich is infiltrated by 

a fluid (labelled h, say). An agent to be sensed is contained witliin tlie fluid. It assumed that the fluid and 
the agent to be sensed have quite different optical properties. Thus, the concentration of the agent within 
the fluid may be gauged by the optical properties of the infiltrated porous material [1, 2, 3]. The optical 
properties used to detect the presence of the agent may be the reflectances or transmittances of the infiltrated 
porous material. Alternatively, if one surface of the porous material were coated with a thin metallic film, 
measurements could be based on the excitation of surface-plasmon-polariton (SPP) waves at the interface 
of the porous material and metal film [4, 5, 6]. 

For example, sculptured thin films (STFs) represent rather promising porous materials for such optical 
sensors [7, 8, 9]. These constitute parallel arrays of nanowires which are grown on substrates by physical 
vapour deposition [10, 11]. By controlled manipulation of the substrate during the deposition process, a 
range of nanowire shapes can be achieved. Thereby, the multiscale porosity of such STFs can be tailored 
to order, to a considerable degree. Additionally, since STFs can fabricated from a wide range of organic 
and inorganic materials, a wide range of optical properties for the porous material can be delivered [12, 13]. 
Chiral STFs are especially interesting for optical sensing applications, as these support the circular Bragg 
phenomenon, courtesy of the helical nature of their nanowires [10]; furthermore, they also support more than 
one mode of SPP wave [14, 15, 16] which may be usefully exploited for sensing [17]. 

In the design of such an optical sensor, what values should one choose for the optical properties of the 
porous material and infiltrating fluid, in order to maximize sensitivity? What value should one choose for 
the porosity, and what shape should one choose for the pores, in order to maximize sensitivity? These 
are the questions that we address here. We do so by considering the simplest scenario wherein the porous 
material and infiltrating fluid are both made from lossless, homogeneous, isotropic dielectric materials, 
characterized by relative permittivities e" and e'', respectively.^ The infiltrated porous material is regarded 
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^Thc entire analysis presented herein may also be applied to lossless, homogeneous, isotropic magnetic materials, by replacing 
relative permittivities by the corresponding relative permeabilities throughout. 
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as a homogeneous composite material (HCM) , which is a reasonable approximation provided that the linear 
dimensions of the pores are much smaller than the wavelength(s) involved. Thus, in the case of optical 
sensors operating at visible wavelengths, we have in mind pore linear dimensions ^ 38 nm for the smallest 
values of e°-'^ considered and ^ 10 nm for the largest values of e"''' considered. The infiltrated porous material 
may be either isotropic or anisotropic depending upon the shape of the pores. We use the well-established 
Bruggeman homogenization formalism to estimate the relative permittivity dyadic of the infiltrated porous 
material, namely e^^ [18, 19]. The Bruggeman formalism has recently been implemented to study the 
prospects of infiltrated STFs as optical sensors [20], based on both changes in reflectance/transmittance [21] 
and SPP wave excitation [17, 22]. 

Regardless of whether the sensor is based on changes in reflectance/transmittance or the excitation of SPP 
waves, the sensitivity of the sensor depends crucially on how much the optical properties of the infiltrated 
porous material change in response to changes in the optical properties of the infiltrating fluid. Thus, the 
derivative de^^ /de^ is a key indicator of sensitivity. In the following we explore how this derivative varies 
as a function of the porosity, the pore shape and the relative permittivities of the infiltrating fluid and the 
porous material. 



2 Homogenization theory 

Within our homogenization framework, the pores are all assumed to have the same shape, which is spheroidal 
in general. These spheroidal pores are randomly distributed but identically oriented. The surface of each 
spheroid relative to its centroid is prescribed by the vector [19] 

Ts{6,cf>)=Vg-H0,^), (1) 

with f being the radial unit vector originating from the spheroid's centroid, specifled by the spherical polar 
coordinates 9 and cj). The linear dimensions of the spheroid, as determined by the parameter r], are assumed 
to be small relative to the electromagnetic wavelength(s). The spheroidal shape is captured by the dyadic 

g=U^l+ (U^^-U^) cc, (2) 

where / is the identity 3x3 dyadic and the unit vector c is parallel to the spheroid's axis of rotational 
symmetry. The linear dimension parallel to c, relative to the equatorial radius of the spheroid, is provided 
by the shape parameter p = U\\/U±. A schematic illustration of such a spheroidal pore is provided in Fig. 1. 

The form of the relative permittivity dyadic e'^'", as estimated using the Bruggeman homogenization 
formalism, mirrors that of the shape dyadic C/. That is, it has the uniaxial form 

|^'- = eri+(ef -ef)cc. (3) 
It emerges as the solution of the dyadic Bruggeman equation [23] 

fa^^+fbi^g, (4) 

where is the null 3x3 dyadic. The scalars fa and ft — I — fa denote the respective volume fractions of the 
porous material and infiltrating fluid. Thus, fi, represents the porosity of the optical sensor. The dyadics 

fi'=(e'l-fO-[l + £-(e'l-|''1]"'' {i = a,b), (5) 

are the polarizability density dyadics of the spheroids in the HCM, while the depolarization dyadic in 
eqn. (5) is given by the double integral [24, 25] 

D^ — [ d<j> [ d9 sine ( 1 | [/"^ • f r • C/"^ (6) 
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It may be expressed in the uniaxial form 

£ = £>_lI+ (D|| -£)_l)cc, 

with components 

Du 



^r||(7), 



1 



r±(7), 



wherein the terms 



are functions of the scalar parameter 
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de 
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The double integrals on the right sides of eqns. (10) and (11) may be evaluated as 
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Notice that the anomalous case 7 < which represents a hyperbolic HCM [26] is excluded from our consid- 
eration. 

The dyadic Bruggeman equation (4) yields the two nonlinear scalar equations 



-fa 



^Br 
t t I 



fa + 



1 + D_L (e« - ef '■) l + D±{eb-ef'- 



.fb = 0, 



-fb = 0, 



(15) 



(16) 



which are coupled via -D±.||- Using standard numerical techniques, tliis pair can be solved for ej^'' and e^^. 

Let us turn to the dyadic derivative which provides a measure of the sensitivity of the porous optical 
sensor under consideration, namely 



de'' de'' 



(17) 
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Before proeeeciing further, we observe that the corresponding derivatives of the depolarization dyadic com- 
ponents may be expressed as 



dD 
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~dJ' 
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and derivatives 
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Next we exploit the scalar Bruggeman equations (15) and (16). Their derivatives with respect to e^ may be 
written as 
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with 



/3n = t.n(e«-ef)(e^-ef)+D||(2.f ^e»-e 

A. = -i2(e'^-er)(e^-ef), 

/321 = r.2i(e"-6f)(e^-ef) , 

/323 = h + D^ (e» - ef ) . 





(28) 



(29) 
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(31) 
(32) 
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Thus, the sought after derivatives of ejj and ej? finally emerge as 



/3ll/322 — /3l2/321 
/32l/?13 — /3ll/323 
/3ll/322 — ,^12/321 



da; 



(34) 



(35) 



3 Numerical investigations 

The consequences of the theory presented in the previous section are illustrated here by means of some 
numerical examples. We begin with the simplest case in §3.1 wherein the pores arc spherical and the 
infiltrated porous material is accordingly considered to be an isotropic HCM. Then the effects of anisotropy 
are considered in §3.2 wherein the pores are taken to be spheroidal in shape. For the purposes of these 
numerical calculations, our attention is restricted to relative permittivity values which, at optical frequencies, 
are attainable either using naturally-occurring materials or currently-available engineered materials. In §4 
the results of implementing relative permittivity values which lie beyond the reach of present-day technology 
are commented upon. 

3.1 Spherical pores 

If the pores are spherical (i.e., p = 1) then the relative permittivity dyadic characterizing the infiltrated 
porous material reduces the the scalar form e^'^ = e^^I_ with e^^ = ej^'' = e^'^. For e" G {1.5,5,15}, 
the Bruggeman estimate e^^ and its derivative de^^/de'' are plotted versus E (lj3) and ft, E (0,1) in 
Fig. 1. We see that when e" = 1.5, the Bruggeman estimate e^^ varies approximately linearly with for all 
fb S (0, 1). However, the relationship between e^^ and e*" becomes increasingly nonlinear as increases. For 
e" = 1.5, the derivative de^"^ /de^ increases in an approximately linear fashion as the porosity ft, increases, 
regardless of the value of e^. However, for = 5 the trend is rather different: here the values of de^^/de'' 
peak at /& ~ 0.7 and the height of this peak rises as decreases. This peak in the value of de^^ /de^ becomes 
more pronounced as the value of e" increases. Indeed, at e" = 15 this peak can be clearly observed even 
when e** = 3. 

3.2 Spheroidal pores 

Let us now explore what happens when the pores are taken to be spheroidal. Accordingly, the HCM 

representing the infiltrated porous material is a uniaxial dielectric material. Following our findings in §3.1, 
we fix e" = 15 in order that the effects of pore shape are more clearly appreciated. The Bruggeman estimates 
e^'jl and their derivatives de^^^/de^ are plotted versus e (1, 3) and e (0, 1) in Fig. 2 for p = 10. Both ej'" 
and fi^'' vary relatively little as increases but both decrease — ej^'" approximately linearly and e^'' more 
nonlinear ly — as increases. The derivative de^'^/de^ increases approximately uniformly as fb increases. 
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for > 1.5. However, the values of de^''' / de^ for < 1.5 are slightly peaked around fb ~ 0.9. The plot 
of de^^/de'^ is similarly peaked, but in this case the peak occurs at fb w 0.5, it is larger in height than the 
de^^/de'' peak, and it extends further into the e'^ > 1.5 region. Also, the height of this de^' /de^ peak is 
substantially larger than the corresponding peak in de^^ /de^ observed in Fig. 1 for e" = 15. 

The pores represented in Fig. 2 are prolate spheroids. The corresponding case of oblate spheroids is 
represented in Fig. 3. The parameters for the plots in Fig. 3 are the same those in as Fig. 2 except that 
p = 0.1. The plot of ej^'' versus and fb in Fig. 3 is very similar to the corresponding plot of e^'" in Fig. 2; 
and likewise for the plots of e^'' in Fig. 3 and ej^'' in Fig. 2. Also, the plots of the derivatives de^'^/de'' 
and def^/de^ versus and in Fig. 3 are similar to the corresponding plots of def^/de^ and de^^ /de^, 
respectively, in Fig. 2, albeit there are qualitative differences in the positions and shapes of the peaks in the 
derivative plots. 

4 Discussion and closing remarks 

An analysis based on the Bruggeman homogenization formalism has provided insights into the sensitivity of 
a generic porous optical sensor. Specifically, for a porous material of relative permittivity e" infiltrated by a 
fluid of relative permittivity e'', we found that; 

• the sensitivity is maximized when there is a large contrast between e'^ and e*; 

• if the contrast between e° and is large, maximum sensitivity is achieved at mid-range values of 

porosity; 

• higher sensitivities may be achieved for close to unity when e" ^ 1, for example; and 

• higher sensitivities may be achieved by incorporating elongated pores. 

In §3 the relative permittivities of the porous material considered were e" e {1.5,5, 15}. These values 
correspond to many common dielectric materials at optical frequencies, with the largest value being close 
to the relative permittivity of silicon, for example. The relative permittivities of the infiltrating fluid were 
taken to be in the range 1 < e*" < 3. This range is physically-realizable, with the largest values corresponding 
to nanocomposite fluids developed for immersion lithography [27] whereas values approaching unity may be 
attained using water vapour [28], for examples. We note that ongoing rapid developments in engineered 
materials are bringing relative permittivity parameter regimes, which were hitherto unattainable, into reach 
[29, 30]. For example, relative permittivities in excess of 50 are now being reported for engineered materials in 
the terahertz frequency regime [31, 32], while engineered materials with positive- valued relative permittivities 
less than unity also appear to be attainable [33, 34, 35]. Accordingly, it is of interest to consider how the 
sensitivities reported here would be effected if rather more exotic parameter regimes were incorporated. In 
further numerical studies (not presented in §3) it was observed that increasing e° beyond 15 results in a 
steady increase in sensitivity. Reducing the value of from unity results in a sharp increase in sensitivity. 
In this context let us make a couple of parenthetical remarks: First, since the entire analysis presented herein 
is isomorphic to the corresponding scenario for magnetic materials (with relative permittivities replaced by 
relative permeabilities throughout), we note that relative permeabilities less than unity can be achieved 
by using diamagnetic materials [36]. Second, the regime e*" < with > (or vice versa) gives rise to 
Bruggeman estimates of the HCM relative permittivity dyadic which are not physically plausible [37] and 
therefore this regime is avoided here. 

The design parameters considered here were the relative permittivities of the porous material and the 
infiltrating fluid, along with the porosity and the shapes of the pores. In the operation of such a generic 
optical sensor — at least in the most straightforward mode of operation — it may be envisaged that the 
relative permittivity of the infiltrating fiuid is a variable quantity (which varies according to concentration 
of the agent to be sensed) whereas the other design parameters remain fixed. Accordingly, the value of the 
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relative permittivity of the fluid should bo carefully chosen such that the sensitivity is maximized over the 
expected range of concentrations of the agent to be sensed. 

While our attention here has been confined to infiltrated porous materials represented as uniaxial di- 
electric HCMs, a straightforward extension of the presented analysis could accommodate biaxial dielectric 
HCMs which represent certain STFs as optical sensors [17, 21, 22]. 

Finally, the study described herein provides a step towards a comprehensive study of porous platforms 
for optical sensing, which incorporates such matters as the absorption/desorption phenomenons that dictate 
the response time and the reversibility of the sensors. 
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Figure 1: Schematic illustration of a spheroidal pore filled with a fluid of relative permittivity e*", embedded 
in a material of relative permittivity e". The spheroid's semi- major and semi- minor axes have lengths [/|| 
and t/_L, respectively, with the semi-major axis being aligned with the direction of c, per eqn. (2). 
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Figure 2: The Bruggeman estimate of the relative permittivity of the infiltrated porous material e plotted 
versus the relative permittivity of the infiltrating fluid e'' G (1, 3) and the porosity ft G (0, 1), for the relative 
permittivity of the porous material e° e {1.5, 5, 15}. Also plotted are the corresponding derivatives de^^ /de'^. 
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Figure 3: The Bruggeman estimate of the relative permittivity parameters of the infihrated porous material 
e^^'ll plotted versus the relative permittivity of the infiltrating fluid e'' € (1,3) and the porosity ft e (0, 1), 
for the pore shape parameter p = 10 and the relative permittivity of the porous material e° = 15. Also 
plotted are the corresponding derivatives de^^.Jde^. 
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Figure 4: As Fig. 3 except that the pore shape parameter p = 0.1. 
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